Current induced magnetization reversal on the surface of a topological insulator 
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We study dynamics of the magnetization coupled to the surface Dirac fermions of a three di- 
mensional topological insulator. By solving the Landau-Lifshitz-Gilbert equation in the presence of 
charge current, we find current induced magnetization dynamics and discuss the possibility of mag- 
netization reversal. The torque from the current injection depends on the transmission probability 
through the ferromagnet and shows nontrivial dependence on the exchange coupling. The mag- 
netization dynamics is a direct manifestation of the inverse spin-galvanic effect and hence another 
ferromagnet is unnecessary to induce spin transfer torque in contrast to the conventional setup. 
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Topological insulator provides a new state of matter 
topologically distinct from the conventional band insu- 
lator—. In particular, edge channels or surface states 
are described by Dirac fermions and protected by the 
band gap in bulk states. Reflecting the topological na- 
ture or the surface Dirac fermion, a number of interesting 
phenomena have been predicted such as the quantized 
magneto-electric effect^i^, giant spin rotation^, magnetic 
properties of the surface state^, magnetization dynam- 
ics^ii, magneto-transport phenomena^^~, and supercon- 
ducting proximity effect with Majorana fermionsi^"— . In 
particular, from the viewpoint of spintronics, the topo- 
logical insulator is an important material to pursue novel 
functionality because spin and momentum are tightly re- 
lated on its surface. In this paper, we utilize this property 
to generate spin transfer torque. 

Spin transfer torque is a fundamental effect in 
spintronics J^^ When a spin current is injected into 
a ferromagnetic layer with a magnetization misaligned 
compared to the polarization of the spin current, the spin 
angular momentum of the injected electron changes upon 
entering and propagating through the ferromagnetic re- 
gion since it follows the direction of the magnetization 
via exchange coupling. This process gives a torque on 
the magnetization, and is called spin transfer torque. 
Consequently, the magnetization can precess^i^ or even 
be switched^i"— . Electrical spin injection has been rou- 
tinely achieved by driving a current through ferromagnet 
with the use of ferromagnet/ferromagnet junction. 

In this paper, we investigate dynamics of the magneti- 
zation coupled to the surface Dirac fermions of a three di- 
mensional topological insulator. By solving the Landau- 
Lifshitz-Gilbert equation in the presence of charge cur- 
rent, we find current induced dynamics of the magnetiza- 
tion and discuss the possibility of magnetization reversal. 
We show that when damping is strong and the current 
induced torque overcomes the anisotropy field, magneti- 
zation reversal occurs. The torque from the current in- 
jection depends on the transmission probability through 
the ferromagnet and shows nontrivial dependence on the 
exchange coupling. The magnetization dynamics is a di- 
rect manifestation of the inverse spin-galvanic effect and 
hence another ferromagnet is unnecessary to induce spin 
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FIG. 1: (Color online) Schematic of the model. By injecting 
charge current on the surface, one can induce magnetization 
dynamics in the ferromagnet deposited on top of a topological 
insulator. 



transfer torque in contrast to the conventional setup. 

We consider magnetization dynamics of the ferromag- 
net deposited on the surface of a topological insulator 
driven by ac charge current fiowing on the surface, as 
shown in Fig. [T] Before proceeding to the explicit cal- 
culation, let us explain the mechanism of the current in- 
duced magnetization dynamics in this setup. Since spin 
and momentum are coupled on the surface of the topo- 
logical insulator, charge current can induce the magne- 
tization, which is the so-called the inverse spin-galvanic 
effect i^ii The emergence of the inverse spin- galvanic ef- 
fect is a direct consequence of the fact that, on the surface 
of topological insulator, the velocity operator is given by 
the Pauli matrices in spin space. The current induced 
magnetization on the surface exerts torque on the mag- 
netization of the ferromagnet when they are noncollinear. 
Therefore, by injecting electric current on the surface of 
the topological insulator, one obtains spin transfer torque 
between the surface Dirac fermion and the magnetiza- 
tion. In sharp contrast, in conventional setup, another 
ferromagnet is necessary to inject spin polarized current 
into the ferromagnet to realize magnetization dynamics. 

Now, let us explain the model to describe the current 
induced magnetization dynamics on the surface of topo- 
logical insulator. The dynamics of the magnetization is 



determined by the Landau-Lifshitz-Gilbert equation: 
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where n is a unit vector pointing in the direction of the 
magnetization, D represents the anisotropy energy with 
easy axis along x direction, ac is the Gilbert damping 
constant, and T is the torque acting on the magnetiza- 
tion. We consider the exchange coupling between the 
magnetization and the surface Dirac fermion of the form 
Hex = m • (T with m = hn and the vector of the Pauli 
matrices in spin space cr. The torque resulting from this 
exchange coupling is then given by 

T = h{cr) xn. (2) 

Since the velocity operator is given by the Pauli matri- 
ces on the surface of a topological insulator, the expec- 
tation value of the Pauli matrices can be represented by 
the charge current which in turn can be obtained through 
the Landauer formula: 



and (uz) = where j, — e, vp, E, and 9 are the cur- 
rent density, the electron charge, the Fermi velocity, the 
n\ Fermi energy, and the angle of incidence, respectively. 
Also, t is transmission coefficient through the ferromag- 
netic region and V represents applied ac voltage given 
by V^ = Vb cos fit. It should be noted that transmission 
coefficient t depends on the magnetization vector m. In 
this way, the problem reduces to the calculation of the 
transmission coefficient. 

To calculate the transmission coefficient, let us explic- 
itly write down the Hamiltonian of the system shown in 

Fig.m 



H — Tivp {kyUx — kx(Jy) + m • <jQ{x)Q{L — x) (5) 
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where Q{x) is the step function and the ferromagnet is 

d9 lil'^sin^ f3) attached to the topological insulator over the region < 

'■ X < L. The magnetization of the ferromagnet is assumed 

to be spatially uniform. The wavefunctions in each region 

d9\t\ cos 9 (4) can be written as 
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where r is the reflection coefficient, E — hvpkp, kx 
kp cos9,ky = kpsm9^ hvpkx = 



E'^ — ml — (hvpky)^, 

and hvpky — fivpky + rux- Due to the translational in- 
variance along the y-axis, the momentum ky is conserved. 
Hence, the common factor e*'^'"*' is omitted above. 

By matching the wavefunctions at the interface x = where 



and L, we obtain the transmission coefficient t: 
4 cos 9hvpkx 
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a = e 



ikpL cos 9 



,ai = e^(^-+™«)i, and a2 = e*(-fe-+™«)'' 



By using the expression of the transmission coefficient. 
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FIG. 2: (Color online) Time evolution of the magnetization 
vector with h/E = 0.1, V = 1, kpL = 100, and aa = 0.1. 
Solid line Q./lof = 0.1. Dotted line Q./ijjf = 10. 



FIG. 3: (Color online) Time evolution of the magnetization 
vector for weaker damping aa = 0.01, h/E = 0.1, kpL — 100, 
and n/ujF = 10. Solid line V = 100. Dotted line V = 1. 



we obtain the torque steming from the current injection 
Eq.(l2]), and finally we solve the Landau-Lifshitz-Gilbert 
equation numerically. In the following, the initial con- 
dition is set to he nx{t = 0) = l,ny{t = 0) = and 
nz{t = 0) = 0. We set uip = D/h and also introduce 



a parameter V 



E^eVo 



2T!hFu 



which measures the applied 



voltage relative to the anisotropy energy since V ex ^. 

In Fig. [21 we show time evolution of the magnetization 
vector for VL/ijjF = 0.1 and 10 with h/E = 0.1, V' = 1, 
kpL = 100, and ac = 0.1. We find the current induced 
magnetization reversal from n^; = 1 to n^j, = — 1. When 
Ux is reversed, the other components Uy and n^ strongly 
oscillate. As 17 increases, the magnetization reversal oc- 
curs at later time. This is because, for large fi, the torque 
from the charge current oscillates very fast compared to 
the characteristic motion of the magnetization. Hence, 
the magnetization would experience an averaged torque 
over many oscillations, which results in small effect due 
to partial cancellation of the net torque. 

Figure [3] shows the time evolution of the magnetiza- 
tion vector for smaller damping constant ciq = 0.01 with 
h/E = 0.1, kpL = 100 and Vt/wp = 10. At i^ = 1, the 
magnetization vector is trapped near the initial condition 
and the dynamics is quite weak (see dotted line). To in- 
duce stronger dynamics, one should increase the torque 
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FIG. 4: (Color online) Time evolution of Ux for aa ~ 0.1, 
V = 1 and Q/iUF = 1 with kpL = 100 and various h/E 
(upper), and with h/E = 0.9 and various kpL (lower). 



arising from the current injection. By increasing the bias 
voltage {V — 100), the dynamics becomes stronger and 
we find strong oscillation of the magnetization, includ- 
ing reversal of the magnetization but the reversal is not 
permanent due to a small damping constant as shown by 
sohd line in Fig. [H 

One may think that for larger exchange coupling, the 
oscillation of the magnetization becomes stronger since 
the torque (Eq.(l2])) seems to become larger. However, 
when the value of rrix becomes large, the Fermi surface 
moves in the fcj,-direction in the ferromagnetic region^. 
Then, the transmission probability becomes very small 
because the number of the evanescent modes increases. 
Therefore, in effect, for large h, the torque does not 
contribute when the magnetization vector points to x- 
direction. Since (cr) lies in the x — y plane, the torque 
(Eq.(l2])) represents an easy axis anisotropy along the axis 
parallel to (cr) in a; — y planeS, but when the magneti- 
zation vector points to a;-direction, the torque becomes 
very small. Thus, for large h, the torque virtually repre- 
sents an easy axis anistropy along y-axis. This indicates 
the deviation of the magnetization from the x-axis for 
large h. We show time evolution of Ux for ac = 0.1, 
V = 1 and Q./ujf = 1 with kpL — 100 and various h/E 
in the upper panel of Fig. 21 and with h/E — 0.9 and 
various kpL in the lower panel of Fig. S) Consistent with 
the above discussion, the magnitude of fix becomes small 
for large ft,, while other components, ny and n^, become 
large. Thus, we find that optimal magnitude of the ex- 
change coupling for the magnetization reversal is in the 
intermediate regime of the exchange coupling (note that 
when ft — > 0, no current induced dynamics occurs since 
T — >■ 0). As L increases, the transmission probability de- 



creases and hence the torque from the current injection 
is suppressed. Thus, the magnitude of the oscillation de- 
creases with increasing L as shown in the lower panel of 
Fig. m However, for smaller ft, since the number of the 
evanescent modes decreases, the effect of the length of 
the ferromagnet L becomes less prominent. 

As for experimental realizability, if we set E =100 meV 
and kp = \ nm"-'^, then the parameter range used in 
this paper is L = 1 ~^ 100 nm and ft = 1 ^ 10 meV. 
Typically, Tiujf = I? ~ 0.1 meV, then V = 1 corresponds 
to eV{) ~ 0.1 meV and the value of O used in this paper 
is in THz regime. These values can be achieved by the 
present experimental techniques. 

In summary, we have studied dynamics of the magneti- 
zation coupled to the surface Dirac fermions of a three di- 
mensional topological insulator. By solving the Landau- 
Lifshitz- Gilbert equation in the presence of charge cur- 
rent, wc have found current induced dynamics of the 
magnetization and discussed the possibility of magneti- 
zation reversal. The torque from the current injection 
depends on the transmission probability through the fer- 
romagnet and shows nontrivial dependence on the ex- 
change coupling. The magnetization dynamics is a di- 
rect manifestation of the inverse spin-galvanic effect and 
hence another ferromagnet is unnecessary to induce spin 
transfer torque in contrast to the conventional setup. 
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